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A generalized Fermi-Bose mapping method is used to determine the exact ground states of six 
models of strongly interacting ultracold gases of two-level atoms in tight waveguides, which are 
generalizations of the Tonks- Girardeau (TG) gas (ID Bose gas with point hard cores) and fermionic 
Tonks- Girardeau (FTG) gas (ID spin-aligned Fermi gas with infinitely strong zero-range attrac- 
tions). Three of these models exhibit a quantum phase transition in the presence of an external 
magnetic field, associated with a cooperative ground state rearrangement wherein Fermi energy 
is traded for internal excitation energy. After investigation of these models in the absence of an 
electromagnetic field, one is generalized to include resonant interactions with a single photon mode, 
leading to a possible thermal phase transition associated with Dicke superradiance. 



(N 



PACS numbers: 03.75.-b,03.75.Mn,42.50.-p 



5-H 

CD 



I 

O 

o 



> 

00 
00 

o 
o 



X 



The rapidly increasing sophistication of experimental 
techniques for probing ultracold gases has caused a shift 
of emphasis in theoretical and experimental work in re- 
cent years, from effective field approaches to more refined 
methods capable of dealing with strong correlations. In 
ultracold gases confined in de Broglie waveguides with 
transverse trapping so tight that the atomic dynam- 
ics is essentially one-dimensional [IJ, with confinement- 
induced resonances P, |2] allowing Feshbach resonance 
tuning Q of the effective ID interactions to very large 
values, such correlations are greatly enhanced. This has 
allowed experimental verification 0, 0, 0] of the fermion- 
ization of bosonic ultracold vapors in such geometries 
predicted by the Fermi-Bose (FB) mapping method [?[■ 
Here the FB mappiiig method used in [7| to solve the 
TG gas and in [3, H, Q to solve the FTG gas is general- 
ized to obtain exact solutions of six models of strongly 
interacting ID ultracold gases of 2-level atoms in tight 
waveguides, with and without resonant coupling to a sin- 
gle mode of the quantized electromagnetic field. 

Consider first some ID gases of ultracold two-level 
atoms with no coupling to the quantized electromagnetic 
field. The procedure for obtaining the exact solutions is 
similar to that used in recent work on two-component 
mixtures of ID gases, but with an important difference: 
In a mixture of dissimilar species there is no symmetry 
requirement on exchange of dissimilar atoms, whereas in 
the present case of a single species with two internal lev- 
els g (ground) and e (excited), the wave function must be 
symmetric (Bose) or antisymmetric (Fermi) under com- 
bined exchange of spatial coordinates and internal states. 
Starting as in [10] from "model states" 4* a/ which are 
ideal Fermi or Bose gases, one can generate the interac- 
tions by a mapping '^m — > ^ = A'^/m where A is ±1 
everywhere and has discontinuities at pair contact which 
introduce zero-range TG and/or FTG interactions. 

Denote spatial coordinates of the N atoms by 
xi, • • • ,xn and their internal state labels by si, • • • , sat 



where each Sj takes on the values g or e. Models I, II, and 
III start from the model states ^^^f °f ideal Fermi gas, 
and models IV, V, and VI from those ^'fjvf of an ideal 
Bose gas: 



*aA/ = (ictfi^-i^<l>u, {xi, se) , = alt"^^i(/)^^ {xe, Si) 

(1) 

where alt denotes an alternant (minus signs in the de- 
terminant det replaced by plus signs). The orbitals 
(/)i,{x,s) are a complete set of energy eigenstates of a 
single-particle Hamiltonian H{x, s) = —^-§^+'v{x, s)-f 
5sg^g + 5se^e appropriate to given boundary conditions 
and external potential v{x,s), satisfying the orthonor- 
mality relation X]s=g e / dx4>l{x, s)4>i,'{x, s) = 5^^' ■ The 
(f)i, in the Slater determinant "^^.^ are a selection of 
N different (f)^, whereas those in the alternant are 
such a selection but with multiple occupancy allowed. 
The terms depending on the internal energy levels Eg 
and eg account for the energy difference ee — eg be- 
tween ground and excited states, "^^i and are eigen- 
states of the interaction-free model Hamiltonian Hm — 
'}2!j=i ^{xj,Sj) with eigenvalue ^v^, and all selec- 

tions of the Vj give all such iV-atom energy eigenstates. 
The Bose ground state is a completely Bose- Einstein con- 
densed state with all N atoms in the lowest orbital v = Q, 
i.e., ^'^jv/ = 4'o{xj, Sj), and the Fermi ground state 

is a Slater determinant of the lowest N orbitals, i.e., 
'^^^lf£=(oi)4'i'ixi, si). Assume harmonic trap- 
ping with the same longitudinal frequency for both in- 
ternal ground and excited atoms, i.e., v{x, s) 
Then if e^g = 



N 
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eg < Nhio, the state with all N atoms 
in internal level g is not the ground state, because the 
energy can be lowered by exciting the atom at the top 
of the Fermi sea to internal level e and moving it to the 
lowest harmonic oscillator level, signalling an instability 
of the putative ground state against internally exciting 
atoms near the top of the Fermi sea and moving them to 
levels near the bottom, which become doubly occupied 
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with one g-atom and one e-atom. If e^g is a hyperfine 
sphtting, this inequahty cannot be satisfied for currently 
achievable values of N , so in the absence of an external 
magnetic field the iV-atom ground states of these models 
consist of atoms in their ground levels g. However, in 
the presence of an external magnetic field in the proper 
direction one has Ceg = ehf ~ ezcoman where enp and 
ezeeman are the hyperfine and Zeeman splittings. This 
can be made arbitrarily small by tuning the magnetic 
field, so a quantum phase transition of this type should 
be experimentally realizable in a TG gas. It also has in- 
trinsic theoretical interest because it is induced by the 
strong interatomic interactions and cannot occur in an 
ideal Fermi gas. 

The energy eigenstates including TG and/or FTG in- 
teratomic interactions are generated by multiplication of 
the ideal Fermi and Bose gas states by appropriate map- 
ping functions A{xi, si, ■ ■ ■ ,xn,sn). The original FB 
mapping solution for the TG gas [7] starts from a model 
state which is an ideal spin-aligned Fermi gas Slater 
determinant, and generates the energy eigenstates of 
a system of bosons with TG (point hard core) interac- 
tions by multiplying 5'^? by a "unit antisymmetric map- 
ping function" 



n ^s^i^j 

i<j<^<w 



xe) 



(2) 



where the sign function sgn(x) is +1 (—1) if x > (x < 
0). This changes the antisymmetry nodes at Xj = xe 
into collision cusps which are the c —>■ +oo limit of the 
collision cusps in the Lieb-Liniger solution for a Bose 
gas with interactions c6{xj — xi) [lH. The FTG gas 
is a "mirror image" of the TG gas, consisting of spin- 
aligned fermions with infinitely-strong zero-range attrac- 
tions which are a zero- width, infinite depth limit of a 
square well of depth Vq and width 2xq, with the limit 
taken such that Vbxg {nh'^)/8Li where /z is the effec- 
tive mass of the colliding pair 0, H, @| . It causes odd- 
wave scattering (ID analog of 3D p-wave scattering) with 
ID scattering length gid = — oo, with the result that 
all energy eigenstates of the FTG gas are obtained 
from corresponding ideal Bose gas states ^ b (the model 
states) by the mapping b ^ f ~ A^b, where A 
is exactly the same mapping ([2]), which now introduces 
sign-changing discontinuities in the FTG states necessary 
to reconcile fermionic antisymmetry with a strong inter- 
action in the zero-range limit Xq ^ Q Inside the 
square well the solution passes smoothly through a zero 
at Xj — = 0, so the discontinuity is an illusion of the 
zero-range limit 0, H, 0, . 

Model I: This is a Bose gas with TG gg and ee in- 
teractions and FTG ge interactions. Its energy eigen- 
states VPf are generated by mapping from the two-level 
ideal Fermi gas model states of Eqs. H]) accord- 

ing to ^'^ = ^^'qtv/ where the mapping function A is 
the simple one ^ of the original TG gas solution f^. 



This generates TG gg and ee interactions, but FTG ge 
interactions. This is easily seen by a simple example. 
Suppose that g and e atoms see the same harmonic trap 
potential v{x, s) = ^miv'^x'^. Then the lowest state "^^m 
with 2 atoms in the ground level g and 1 atom in the 
excited level e is a Slater determinant ^ constructed 
from the three harmonic oscillator orbitals 4'u{x, s) where 
(t)o{x,s) = uo{x)58g, (t)i{x,s) = ua{x)Hi{Q)5s g, and 
02(a;, s) = uo{x)5se^ where Q = x/xosc, Xosc = ^Jh/muj, 
uq{x) = const. e"*^ and _ff„ are Hermite polynomi- 
als, as in ^is'l. Dropping a normalization constant and 
multiplying by A of Eq. ([2]) to obtain one finds 

= Uo[xi)U(i{x2)uo{xz) 

x[-|xi - a:2|sgn(a;i - a::3)sgn(a;2 - X3,)5s^g5s2g5s:^e 
+ \xi - a::3|sgn(a;i - a;2)sgn(x2 - X:i,)5sig5s3g5s2e 
~\x2 - a;3|sgn(a::i - a;2)sgn(a;i - X:i)5s2g5s3g5ste] ■ (3) 

If atoms 1 and 2 are both in level g, then has a TG 
\xi — X2I cusp at Xi — X2, but if atom 1 is in level g 
and atom 2 in level e or vice versa, then has a FTG 
discontinuity at Xi — X2- Similar remarks apply to the 
atom pairs (1,3) and (2,3). This generalizes to eigenstates 
with arbitrary numbers of atoms in levels g and e, 
showing that two atoms both in level g or both in level 
e have TG cusps at their collision points Xj = xi, but if 
one is in level g and the other in level e, there is a FTG 
discontinuity instead. This model is equivalent to that 
of Mousavi et al. 14|, although the equivalence is not 



immediately apparent since the notation in 1J| is very 
different from that used here, and the interactions in 
are defined differently, as TG in the three triplet gg, ee, 
and ge channels and FTG in the singlet ge channel [l5|. 
To see the equivalence note that ([3]) can be rewritten in 
terms of triplet and singlet ge channels as 

= ^Uo{xi)uo{x2)uo{x3){\xi - 0:21 

xsgn(a;i - a;3)sgn(a;2 - X3)6s3g{Ss,g6s2e + Ss2gSsie) 
-\xi - a:3|sgn(xi - X2)sgn(x2 - X3) 

+ \x2 - X3|sgn(xi - a::2)sgn(a;i - 2:3) 

-|-sgn(a;i - a;2)sgn(a;i - a;3)sgn(a;2 - X3) 

X [{Xi - 2X3 + X2)Ss3g{SsigSs2e - Ss2gSsie) 

+ {X3 - 2X2 + Xl)Ss2giSs3gSsie - SsigSgse) 

+ {X2 - 2xi + X3)Ssig{Ss2gSs3e " (Jsgg^sje)]} . (4) 

It follows that the interaction in the triplet ge channels 
is of TG form 17] whereas that in the singlet ge chan- 
nels is of FTG form, but the equal expression ([3]) is much 
simpler due to cancellations in ([4]), does not require sep- 
aration into singlet and triplet channels, yet exhibits the 
TG interaction in the triplet gg channel, which is not 
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evident in the form ([4]). The model represents a two- 
level hybrid TG-FTG gas with infinitely strong repulsive 
gg and ee TG interactions and infinitely strong attrac- 
tive ge FTG interactions, although it was called simply 
a two- level TG gas in [ijj. Its experimental realization 
would be difficult since generation of the FTG interaction 
would require a p-wave ge resonance, and one would have 
to simultaneously create TG gg and ee interactions, al- 
though the latter could perhaps be done by the methods 

of 

Model II: This is a Bose gas with TG gg and ee inter- 
actions but no ge interactions. It starts from the same 
ideal Fermi gas model states '^^j^j of Eqs. ([1]) and gener- 
ates the states ^I'^ with TG gg and ee interactions by a 
mapping = A^^jyj, but the mapping is now a more 
complicated one depending on both spatial and internal 
variables, which is everywhere ±1 and antisymmetric un- 
der exchanges (xj, Sj) ^ {xg, sg): 

A{xi,si;- ■ ■ ;xn,sn) = Yl a{xj,Sj;xe,se) 

a{xj,Sj;xg,se) 

= Ss^,siSgn{Xj - Xt) A- Ss-gSs.e - Ss-eSsig ■ (5) 

This is the same mapping used previously in a model 
of a TG Mach-Zehnder interferometer, where the state 
variable s distinguished between the two interferometer 
arms ■ It generates a TG interaction in the gg and ee 
channels Sg^^sn but no ge interaction. It should be easier 
to realize experimentally than model I. 

Model III: This is a Fermi gas with FTG ge interac- 
tions but no gg or ee interactions. It again starts with 
the ideal Fermi gas model states "^^^m O' but 

uses a different mapping A, which is everywhere ±1 and 
is now symmetric under exchanges (xj^Sj) <-> {xt,se), to 
generate the states Vf^ with FTG interactions: 

A{xi,si;- ■ ■ ;xn,sn) = Y\. a{xj,Sj;x£,Si) 
i<j<e<N 

a{xj,Sj;xe,Si) = Ss^^si 

+ iSsjgSsee - Ss-eSsi,g)sgn{xj - Xl) . (6) 

It generates no interaction in the gg and ee channels 
Ssj,st, but FTG ge interactions. 

Models IV, V, and VI: These all start from the ideal 
two- level Bose gas model states 4*^^/ of Eqs. ([T]). Model 
IV uses the simple mapping ^ and has FTG gg, ee, 
and ge interactions, model V uses the mapping (O and 
has FTG gg and ee interactions, and model VI uses the 
mapping © and has FTG ge interactions. Model VI 
should be the easiest to realize since a p-wave Feshbach 
resonance is required only in the ge channel. 

Interaction with a coherent photon mode: Some thirty 
years ago I suggested the possibility of a cooperative cou- 
pling between Bose-Einstein condensation (BEG) in liq- 
uid ^He and BEG of van der Waals virtual photons into 



a superradiant mode, and constructed a simple model 
exhibiting such a quantum phase transition [l9|. Such 
equilibrium superradiance has never been observed (nor, 
I believe, searched for) in superfluid ^He, but in recent 
years coupling of matter waves and a superradiant pho- 
ton mode have been observed in ultracold Bose gases 

20l . 21 1 , and it has been pointed out that essentially the 
same phenomena can also occur in ultracold Fermi gases 

S 23I. In fact, it had been shown long ago by Hepp and 
Lieb [2J] that BEG of the atoms is not necessary for co- 
operative coupling of internal and translational degrees 
of freedom in a multilevel gas via a superradiant photon 
mode. Motivated by these results, model I above will be 
generalized by adding resonant coupling to a single mode 
of the quantized electromagnetic field [25| . 

Start with the atom-field interaction Hamiltonian in 
the electric field gauge and rotating wave approximation: 



H 



atom-field 



fcA 



(deg • ekxSfbkxe 



H.c). Here bk\ and are annihilation and creation op- 
erators for photons with wave vector k and polarization 
A, Sj' and S~ ~ (S'j')^are raising and lowering operators 
for the internal levels in the usual spin ^ representation, 
deg is the transition dipolc moment from level g to level 
e, e^A are unit polarization vectors, and the allowed k 
are integer multiples of 27t/L where L is the length of a 
ID microwave cavity resonant with the hyperfine transi- 
tion g e. Retaining coupling only to a single mode of 
wavelength A = ^ resonant with the hyperfine transi- 
tion, one obtains a Hamiltonian generalizing that of mod- 



els I-III: H 

1 



'2m aJf 



-v{xj) 



H.C.)] where q = 



Bq, a constant term 



N„ 

{hcq + Eg 



= 

' ee) has 



been dropped, and the external potential is taken here to 
be that of a longitudinal harmonic trap, v{x) = ^mcj^x'^. 
To obtain the ground state with both coupling to the 
electromagnetic field and TG and/or FTG interatomic 
interactions, one can first let the above Hamiltonian act 
on the space of ideal Fermi gas model states ^'f^/, find 
the ground state ^'^^^ in that space, and then find the in- 
teracting ground state by application of the mapping 



The length of typical cigar-shaped ultacold gas traps 
of high aspect ratio is many orders of magnitude smaller 
than relevant microwave wavelengths, so it is an excellent 
approximation to replace e"^^^ by unity (dipole approxi- 
mation for the whole A^-atom system) , leading to a zero- 
order Hamiltonian Hq = i?trans + ^^Dickc where the trans- 
lational Hamiltonian is i/trans = EiLiI"^^ + ^i^j)] 

and the remainder is a resonant Dicke model, -ffDickc = 
eeg{Ng + S') + ilL-^'^Y.AS+bq - H.c.)] where = 



Sj and the resonance condition hcq — e^g has been 
inserted. There is strong statistical field-atom coupling 
since the eigenstates of //trans are those of an ideal Fermi 
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gas which must be antisymmetric under combined space- my research, 
spin exchange {xj, Sj) (xi, Si). 

The ground state of i?Dicko in the thermodynamic hmit 
is the direct product of a field-independent state and 
a Glauber coherent state for the electromagnetic field 
2^ 27 1 . This amounts to replacement of the photon anni- 
hilation and creation operators by c-numbers Pg and /?* . 
The ground state energy is invariant under a gauge trans- 



formation of the first kind (]„ — > /3„e* with 9 arbitrary; a 



convenient choice here is (3q = —iVNfq where fq is real, 
nonnegative, and independent of N. The corresponding 



reduced Hamiltonian is Ho = e^gN fq + X^^Li l 



2m dx'-^ 



vixj) + eegS^ + i^qfq{S+ + SJ)] whcrc pq = N/L. 
The single-atom energies and eigenstates </)^ can be 
found by a 2 X 2 matrix diagonalization for each j , with 

the results = nfiw ± \ ^^ig + ^I'^Pqfq and 0^(a;, s) = 

Un{x)w^{s) where n — j~l = 0, 1,2,---, Un{x) 
the harmonic oscillator eigenfunctions, and 

= -i) = (6eg ± ^elg+4j^pqfq)/2j./^q. 

The quantum phase transition previously described in 
the presence of an external magnetic field B will still oc- 
cur, but now at a higher value of B since — > Ceg- 
The above assumes that the trapped ultracold gas 
is bathed in a constant microwave field which is sup- 
plied and controlled externally. If it is instead con- 
tained in a perfectly refiecting cavity with no exter- 
nal microwave source, one can investigate the possibil- 
ity of a Dicke thermal phase transition to a state with a 
self-generated superradiant microwave field. Minimizing 
the grand canonical free energy —P~^\nTi e~^''^°^'^^^ 
of the ideal Fermi gas of models I-III with respect to 
[27! yields the following condition for determin- 



ing 



2Neegfq = 



7 Pq 



' - S-o(5r7 - 9i) Where 

= [1 + e''^'" ^'^^]^^ and p, = ep determined from 
^^=o(.9n + 5n ) — ^- Even for the much simpler case 
of fixed atoms treated classically (no Fermi sea), solu- 
tion for the superradiant phase transition is nontrivial, 
so I shall not proceed further now. However, by compar- 
ison with that case [2^, |23| it is reasonable to conjecture 
that there is no phase transition for small values of the 
ratio 7y^/eeg but a thermal transition to a superradi- 
ant phase for large values of this ratio. Recalling that Cgg 
can be made arbitrarily small by tuning an external mag- 
netic field B, it seems likely that this model will exhibit 
a superradiant phase transition for sufficiently large B, 
which may be coupled to the previously described coop- 
erative rearrangement of the translational ground state. 
Note that arguments against existence of a superradiant 
transition in real atomic systems J2§\ do not apply here, 
since Ceg includes a negative Zeeman shift due to B. 
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